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Abstract:This paper presents some results on some aspects of the two-level 
atom interacting with a single-mode with the privileged field mode being in the 
squeezed displaced Fock state (SDFS). The exact results are employed to perform 
a careful investigation of the temporal evolution of the atomic inversion, entropy 
and phase distribution. It is shown that the interference between component states 
leads to non-classical oscillations in the photon number distribution. At mid revival 
time the field is almost in the pure state. We have briefiy discussed the evolution 
of the Q function of the cavity field. The connection between the field entropy and 
the collapses and revivals of the atomic inversion has been established. We find that 
the phase probability distribution of the field refiect the collapses and revivals of the 
level occupation probabilities in most situations. The interaction brings about the 
symmetrical splitting of the phase probability distribution. The general conclusions 
reached are illustrated by numerical results. 
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1 Introduction 



The entropy of a radiation field is one of the canonical problems of statistical physics 
and has attracted much attention in the past. In recent years much attention has 
been focused on the properties of the entanglement between the field and the atom 
and in particular the entropy of the system [1-11]. The authors in [2-4] have shown 
that entropy is a very useful operational measure of the purity of the quantum 
state, which automatically includes all moments of the density operator. The time 
evolution of the field (atomic) entropy refiects the time evolution of the degree of 
entanglement between the atom and the field. The higher the entropy, the greater 
the entanglement. 

The concept of the photon in the quantum theory of a radiation field has been 
built on the number (Fock) state |n). However, another important state is the co- 
herent state which is a linear superposition of all \n) states with coefficients chosen 
such that the photon number distribution is Poissonian. It may be defined by the 
action of a displacement operator D{q) on the vacuum state. This state has been 
extensitevly studied [12-14]. On the other hand the squeezed state is one of the 
non-classical states of the electromagnetic field in which certain observables exhibit 
fiuctuations less than in the vacuum state [15-17]. It is defined by the action of the 
squeeze operator S{z) on the coherent state [15-17]. Squeezed displaced Fock states 
(SDFS) have been introduced and difi^crcnt aspects of these states such as squeez- 
ing and photon statistics have been investigated [18-25]. These states generalize 
two- photon coherent states [17] (squeezed coherent states), squeezed number states 
[18], and displaced Fock states [26-28]. They exhibit both number and quadrature 
squeezing. Recently the creation of nonclassical states of motion of a trapped ion 
such as Fock states, coherent states, squeezed states and Schrodinger-cat states have 
been reported [29-31]. In these experiments an ion is laser-cooled in a Paul trap to 
the ground harmonic state. Then the atom is put into various quantum states of 
motion by applications of optical and electric fields. That moved the study of these 
states from the academic realm to the world of experimentation. This motivated us 
to study the interaction of these states with a two-level atom. 

A stochastic formulation of quantum mechanics involves, basically, two inter- 
related problems. These are the determination of the probability functions of the 
density operator, p, and the establishment of the proper correspondence between 
quantum-mechanical observable, and ordinary functions in phase-space. Attempts 
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in this direction ran into the difficulty of deahng with quasidistributions, [32-33]. 
The three types of quasiprobabihty distribution functions P, Q and W (for normal, 
antinormal and symmetric ordering respectively) are very important in quantum 
optics [32-33]. 

Recently Pegg and Barnett [34, 35] have introduced a new hermitian phase for- 
malism which successfully overcomes the troubles inherent in the Susskind-Glogower 
phase formalism and enables one to study finer detials of the phase properties of 
quantum fields. Such quantities as expectation values and variances of the hermitian 
phase operators or phase distribution functions arc now available for investigation 
[36, 37]. One of our interests is to investigate the phase properties here. 

The material of this paper is arranged as follows: In section 2, we review a few 
concepts of squeezed displaced Fock states (SDFS's). In section 3 we introduce the 
model and write the expressions for the final state vector at any time t > 0. We 
discuss the field entropy in section 3.1. By numerical computations, we examine 
the influence of the SDFS's on the field entropy evolution and entanglement of the 
atom and the field. We analyse oscillations in the photon number distribution of 
the cavity field in section 3.2. Section 3.3 discuss the phase probability distribution 
through the framework of Pegg-Barnett's definition of the Hermitian phase operator. 
We present the evolution of the Q function for the SDFS's in section 3.4. Finally, 
summary and remarks are presented. 



2 Squeezed displaced Fock states (SDFS's) 

The SDFS's are generated from the number state |m) as shown below. These states 
have been studied extensively in literature, because of their interesting nonclassi- 
cal properties and expected prospective applications in optical communication and 
interferometery [14]. The SDFS is defined by [18-25] 

\ao, Zjin) — D{aQ)S{z)\'m) (1) 

where the displacement operator D{ao) ,( with ao a complex parameter that repre- 
sents the magnitude and angle of the displacement) [12-13], and the squeeze operator 
S{z) are given by [14-17] 

D{ao) = exp(Q;oa^ — ao*a), 
S{z) = exp[|a2-|a+% (2) 
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where z — re*"^ and r is known as the squeeze parameter and indecates the direction 
of the squeezeing, with a (a+) are the annihilation (creation) operators of the field. 

It is easy to calculate the average value of the number operator, a+a, in the state 
|Q;o,^,m), by using the above relations, thus 

(a+a) = + li/^m + + |Q;o|^ (3) 

where = coshr, and u = exp{i(f)) sinhr. 

The photon number distributions P„ for |ao,2;,'n2) is equal to the square of the 
absolute value of the matrix elements {n\a, z,7n). The analytical expression for 
{n\a,z,m) is given by [19], 

(n| ao,z,m) = {^y/'ij-.r^' exp{-\-^ + f^aj) 



^Z^i=0 \i ) {n-i)\ Lv / 



X ] -f^n-i [ (2,/^) 1/2 ] Hm-i [ (_2z/* Jl) 1/2 ] > (4) 

where cio = /^Q^o + ^^^o*) and Hn{x) stands for the Hermite function of order n [38]. 
Therefore, the photon number distribution P„ is given by 

Pn = \ {n\ao,z,m)f, (5) 

The scalar product of two different SDFS states (cti, zi, mi|Q;25 ^2, ''^2) is very 
useful in the representation with SDFS basis. We can obtain it with the help of the 
completeness of coherent states as 

(ai,2;i,mi|a2,2:2,r"2) = \ j {oii, Zi,mi\(3){P\a2, Z2,m2) d'^(3 

By using equation (4), 



( ai,zi,mi\a2,Z2,m2)^ 
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X 



mi —2s' 



-2//ii/i)V2j (^2-2s-r)!(mi -2s-r)! 



Z^2 



mi — 2s— r 



mi— 2s— r 



X 



(q;o2 - aoi) - ^(ttoi - "02) 



rl"5l -"02, 
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m2— 2s— r 
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- ^("02 - tt0l)l 



where = /U2/ii — u 2^1/ 1^21^1 ^-iid aoi = j^iai — z^io;*, ao2 = Ai2'^2 — i'2(y.2- An 
equivalent result (however not having the Hermite polynomials in the summand) 
has been given recently [24]. When mi = 1712 = the above equation gives the 
result in [16], for the squeezed state. But when = 1/2 = and /ii = = ^ then 
this equation becomes 



mm(mi,m2) 



{ai,mi\a2,m2) = 



E 



X- 



m2\m\\ 



(ai|a2) 



r!(m2 — r)!(mi — r)\ 

which is the result of, [28,18], for the displaced Fock state. Other special cases follow 
in a straightforward way. 



3 The quantum dynamics 

We consider the Hamiltonian for the one-photon Jaynes-Cummings model. It de- 
scribes the interaction of a single-mode quantized field with a two-level atom via a 
one-photon process. The Hamiltonian of the system in the rotating-wave approxi- 
mation is written as 



if = -a;oCr2 -|- a;a^a -|- A(a(7+ -|- aV_), (6) 

where (a) is the creation (annihilation) operators for the photon of frequency 
tu, and A describes the coupling to the atomic system. The two level atom with 
transition frequency lOo is described by the Pauli raising and ( lowering ) operators 
(T+, ((T_) and the inversion operator cr^. with the detuning parameter A = a; — lOo- 
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The initial state of the total atom-field system can be written as 



I V(0) >= V'(O) >/ ® I V'(O) >a= J] I e >, (7) 

n=0 

means that the atom starts in its excited state, the field is assumed to be initially 
in the squeezed displaced Fock states(SDFS) Qn —< n \ ao,z,m > given by Eq. 
(4). The solution of the Schrodinger equation in the interaction picture i.e the wave 
function of the system at any time t > is given by 

oo 

I >= Yl (^"(^) I ^' ^ > I n + 1, ^ >) , (8) 

n=0 

where the coefficients and S„ are given by the formulae 

/ iA sin Xti/nX 

A-n[t) = Qn [cos XtUn - — j, (9) 



„ / V / sin XtUn , V 

Bn{t) = -iqn^n + 1 ^, (10) 



A2 



^n=\/7T2+^+l (11) 



With the wave function | V'(^)) found, any property related to the atom or the field 
can be calculated. The reduced density operator of the field of the system can be 
written as Pf{t) = Tratom\ i^{t)){ip{t) |, 

oo 
n,m=0 

|n + l)(m+l I (12) 

= \C{t)){C{t)\ + \S{t)){S{t)l 

i A sin XtVr, 



^W) = > ^ qn ( cos XtVn " ^ j 



n) 



n=0 

oo 



S{t)) = _i^5„V^^^^^|n + l). (13) 

n=0 
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3.1 Atomic Inversion 



Using equation (8) we can evaluate the time evolution of the atomic inversion W{t) — 
{iP{t) I a, I iP{t)) 

oo 

mt) = E {MtWnit) - Bnmiit)) , (14) 

n=0 

with Pn = l^np of equation (5) describing the distribution of the photons in the 
state SDFS. We start our analysis by taking the squeeze parameter r — 1 and for 
different values of m (m = 0, 1, 2). Since the resulting series cannot be analytically 
summed in a closed form, we will evaluate them numerically. In figure la we plot 
the atomic inversion as a function of scaled time Xt for the squeeze parameter r — 1 
and with zero value of {m — 0). 

We know that the collapses are caused by the dephasing of the various terms in 
the sums in equation (14). Thus we can calculate the time in which the revivals 
will occur by estimating the time that neighbor terms in the sums will be in phase 
again (for n ~ n): Tr(2A-\/^ + 1 — 2A-\/^) ~ 27r. This argument is true in the case 
of coherent state and can be applied for the squeezed coherent state (i.e, m — 0), 
as can be seen in figure la, but it cannot be applied for the SDFS. For nonzero m 
(m = 1) figure lb, not only the amplitude of Rabi oscillations, but also the time 
average of the inversion is affected. In figure Ic we have increased m to 2, and taken 
all other parameter as in figure la. One observes that for longer time the inversion 
shows small oscillations around zero but in quite irregular manner. 

3.2 Entropy of the cavity field 

Employing the reduced field density operator given by Eq. (12), we investigate the 
properties of the entropy. The quantum dynamics described by the Hamiltonian (6) 
leads to an entanglement between the field and the atom. We use the field entropy 
as a measurement of the degree of entanglement between the field and the atom of 
the system under consideration. In order to derive a calculation formalism of the 
field entropy, we must obtain the eigenvalues and eigenstates of the reduced field 
density operator given by Eq. (12). A general method to calculate the various field 
eigenstates in a simple way can be found in [4] . By using this method we obtain the 
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Figure 1: The atomic inversion W{t) for the output field against the scaled time 
Xt. The values of the parameters are assumed as follows, the squeeze parameter 
r — 1, the initial position (3,0), i.e ckq = 3, and the direction of squeeze is along 
the x-direction, i.e (p — 0. While the number of photons for the input state m takes 
different values, where (a) m — 0, (h) m — 1, and (c) m — 2 . 
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eigenvalues and eigenstates of the reduced density operator, 

Xf{t) = {C{t) I C{t)) ± bMtO] I {C{t) I S{t)) I (15) 

= {S{t) I S{t)) ± exp[±^] I {C{t) I S{t)) I, (16) 

I V'/ (^)) = -^=^^={exp[{z<P ± ^)/2] I C{t)) (17) 
j2A^(t) cosh(^) 

±exp[-(z0±^)/2] |5(t))}, (18) 



where 



{C{t) 1 


Cit)) - 


(5W 1 


S{t)) 


2 


1 (c(t) 


s{t)) 





g^sinh-^ P"^''^', ;;vri o7^ ^'M - (19) 



We can express the field entropy Sf{t) in terms of the eigenvalue X^{t) of the reduced 
field density operator, 

Sf{t) = -[X-^{t)lnX-^{t) + X]{t)lnX]{t)]. (20) 

It does not appear possible to express the sums in equation (20) in closed form, 
but for not too large n, direct numerical evaluation can be performed. On the basis 
of the analytical solution presented in the previous section, we shall examine the 
temporal evolution of the field entropy. It should be emphasized that in computing 
all infinite series for the atomic wave function \ip{t)),we have invoked mathematically 
sound truncation criteria. To ensure an excellent accuracy the behavior of the 
field entropy function Sf{t) has been determined with great precision. For regions 
exhibiting strong fluctuation a resolution of 10'^ point per unit of time has been 
employed. The time t has been scaled; one unit of time is given by the inverse of 
the coupling constant A. 

Wc display the evolution of the flcld entropy for different values of mean photon 
number m and the squeeze parameter r is taken to equal (r = 1). In our compu- 
tations, we have taken the displacement parameter cto = 3, A = 0. In the case of 
an initial field with r = (coherent state), we already know [39-40] that the field 
evolves and (nearly) returns to a pure state only at half of the revival time. As fur- 
ther analysis showed, at that time the field is indeed in a superposition of coherent 
states. However, if we initially prepare the field in a squeezed state, as we see in 
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Figure 2: Entropy Sf{t) of the output field against the scaled time Xt. The values 
of the parameters are assumed as follows, the squeeze parameter r — 1, the initial 
position (3,0), i.e ckq = 3, and the direction of squeeze is along the x-direction, i.e 
(j) — 0. While the number of photons for the input state m takes different values, 
where (a) m = 0, (b) m = 1, and (c) m — 2 . 
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figure 2, (r = 1) that its entropy reaches a minimum at approximately the revival 
time corresponding to an initial coherent state with amplitude a. Note that the 
revival time in the case of an initial squeezed state is \t = 2ti{\ a p + sinh^ r)^/^. As 
we see in figure 2a, in the squeezed state case, the value of the entropy is almost the 
same at the revival time and at half of the revival time, and because they correspond 
to minima, the states of the field are less mixed at these times. However, near the 
minima the behavior of the entropy is different for the two cases. While the min- 
imum is reached in a smooth way for the coherent state, it becomes oscillatory for 
the squeezed. The programme described near Sf(t) has been carried out for several 
parameter sets, including those covered by figure 2. It may seem rather surprising 
to have a pure state at a time different from half of the revival time, but this is of 
course due to the nature of the squeezed states. Also, by increasing the squeezing 
parameter r, the field entropy becomes increasingly irregular and with values charac- 
teristic of a mixed state, Spit) — 0.7. This is in qualitative agreement with the fact 
that the atomic response for the field initially prepared in a squeezed vacuum state 
is very similar to when it is prepared in a thermal (mixed) state [41]. To visualize 
the influence of SDFS in the held entropy we set different values of m (m = 1, 2) 
and the squeeze parameter r = 1, and all the other parameters are the same as in 
flgure 2a. The outcome is presented in figure 2b, c. One can distinguish between 
two stages of evolution, each of which has been pictured separately. We see that 
the amplitude of the oscillations as well as the revival time becomes smaller and we 
have more revivals in the same time, the collapse time decreases. It is evident that 
the field and the atom are in the entangled state when m increases further. In next 
section we turn our attention to interesting non-classical phenomenon emerging as 
a direct consequence of quantum interference between component states of the field. 
Namely, we will analyse oscillations in the photon number distribution of the cavity 
field at different values of m. 

3.3 Photon number distribution 

In direct detection one counts the number of photons in the field mode of interest. 
The probabilty for finding n photons: at time t > 0, is given by 

P{n,t) = \{n\m)? = {n\pf{t)\n) (21) 

from which we can easily find the photon number distribution of the cavity 
field in the one-photon JCM. As the cavity field starts to interact with the atom 
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Figure 3: Plots for the photon number distribution P{l,t) of the initial field. The 
values of the parameters are assumed as follows, the squeeze parameter r — 1, the 
initial position (0.5,0), i.e ckq = 0.5, and the direction of squeeze is along the x- 
direction, i.e (f) — 0. While the number of photons for the input state m takes 
different values, where (a) m = 0, (b) m = 1, and (c) m — 2 . 
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the initial photon number distribution P„ =| ^ starts to change. Due to the 
quantum interference between component states the oscillations in the cavity field 
become to be composed of two component states. Even though the field entropy 
is not equal to zero these two component states partially interfere which results in 
some oscillations of the photon number distribution. In figure 3 we plot the photon- 
number distribution of the cavity field in the case of nonzero squeezing parameter r = 
1 and for different values of m. We note that the amplitude of the photon-number 
distribution decreases as m increases see figure 3c (where we have set m=2). For m = 
0( figure 3a), the photon number distribution resembles a Poissonian distribution. In 
the large values of m, the amplitude of the oscillations is affected. As compared with 
the case m = 0, the locations of the maxima have moved to the right as m increased 
further. For the entropy the characteristic sequence of minima is there again (figure 



3.4 Phase distribution 

Recently, Barnett and Pegg defined a Hermitian phase operator in a finite dimen- 
sional state space [34-36]. They used the fact that, in this state space, one can 
define phase states rigorously. The phase operator is then defined as the projection 
operator on the particular phase state multiplied by the corresponding value of the 
phase. The main idea of the Pegg-Barnett formalism consists in evaluation of all 
expectation values of physical variables in a finite dimensional Hilbert space. These 
give real numbers which depend parametrically on the dimension of the Hilbert 
space. Because a complete description of the harmonic oscillator involves an infinite 
number of states to be taken, a limit is taken only after the physical results are 
evaluated. This leads to proper limit which correspond to the results obtainable in 
ordinary quantum mechanics. It can be used for investigation of the phase properties 
of quantum states of the single mode of the electromagnetic field [34-37]. 

The Pegg-Barnett phase distribution P{r}) is defined through an infinite sum 



2c). 



[34,35]: 



(22) 



1,3=0 



where the density matrix given by 



(23) 
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Figure 4: Shows the plot of phase distribution P{r),t) as a function of the scaled 
time Xt and rj. The values of the parameters are assumed as follows, the squeeze 
parameter r — 1, the initial position (3,0), i.e ckq = 3, and the direction of squeeze 
is along the x-direction, i.e = 0. While the number of photons for the input state 
m takes different values, where (a) m = 0, (b) m = 1, and (c) m = 2. 

Figure. (4) shows the time evolution of the phase probability distribution P{ri,t) 
for r = 1 and for various values of m. When m equal to zero, it is remarked that 
P{rj,t) exhibits symmetric splitting as At varies as shown in figure (4a) This is 
the counterrotating behaviour observed earlier [37]. When Xt — 0, P{r],t) has a 
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single-peak structure corresponding to the initial coherent state. The peaks are 
symmetric about 77 = so that the mean phase always remains equal to zero. The 
time behaviour of the phase probability distribution carries some information about 

the collapse and revival of Rabi oscillations [39]. When the phase peaks are well 
separated the Rabi oscillations collapse and each time as the peaks meet { at 1] = 
and/or ±7r ) they produce a revival see figure 4a. It is further noted that the height 
of the peak change as time develops in contrast to the case of the coherent input 
[39]. When m 7^ 0, the situation is completely changed, as we observe from figure 
(4b-c). It is seen that when m = 1, thee are two peaks, one with small amplitude 
compared with the other peak, whose rate of shift becomes faster when plotted in a 
phase space as in [35] . It is observed that the symmetry shown in the the case when 
m = for the phase distribution is no longer present once the new state is added. 
The peaks are split but the two split peaks move with different rates. The one with 
the slower rate faces damping while the faster peak changes in the amplitude as time 
develops. 



3.5 Evolution of the Q function 

In the previous section we discussed a particular aspect of the atomic dynamics 
(collapses and revivals) in the one-photon model with the initial field prepared in 
the squeezed displaced Fock states. Now we are going to try to understand better 
the behaviour of the system by focusing our attention on the field dynamics by 
studing the quasi- distribution function. The first step to be taken is the calculation 
of the reduced density operator of the field Pf{t) (see equation (12)), we get for the 
Q function 

Q{a,t) = ^{a\pf{t)\a) (25) 

where \a) is a coherent state. More than just a theoretical curiosity, Q{a,t) can be 
detected in homodyne experiments [33] . It has the form 

^/ X 1 / , ,2x / ,(a*)"(a)"* , , 

Q{a,t) = -exp(-|an Pumit V l-^ (26) 

The Q function is not only a convenient tool to calculate expectation values of anti- 
normally ordered products of operators, but also gives us a new insight into the 
mechanism of interaction in the model under consideration. 
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Figure 5: The Q function for the output of the JCM driven by the SDFS. The 
parameters are assumed as follows: The squeeze parameter r = 1, = 0, «o = 3, 
A = 0, and m = 1. While the interaction time takes the values: (a) t = 0] (b) 
t = tji/2; (c) t = tR. Here X = Re{a) and Y = Im{a). 
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In figure 5 we have sketched the Q-function distribution functions for the output 
field states. We here show this for the squeeze parameter r — 1, ag = 3, m — 1 and 
for different values of the interaction time. We see that the Q function corresponding 
to the initial squeezed coherent state, [16], figure 5a bifurcates due to the quantum 
interaction between the field and the atom. At At = the Q function has a single- 
peak structure, but with At = Tr/2 we see that the Q function is composed of two 
well separated components, mean while if At = Tr, three components have been 
seen figure 5c. We show that at t = the peak of SDFS is observed [19]. For t > 
the two peaks split into two sets of counter-rotating peaks during the collapse. At 
longer times the Q-function is spread out over an angular region in the xy-plane 
as shown in figure 5. If we combine this observation with the fact that the field 
entropy at this moment is almost equal to zero we can conclude that the cavity field 
at m = is in a pure state. On the other hand at m = 1 the entropy reaches its 
maximum and so in spite of fact that the Q function is composed of two parts the 
cavity field is in a statistical mixture state. 



4 Summary and concluding remarks 

The well-known Jaynes-Cummings Model gives an exactly solvable model of a two- 
level atomic system in interaction with a radiation field. We considered the inter- 
action with the field initially is the SDFS. We have been studied in this work the 
three aspects, the first is the dynamics of the atomic inversion, the second is the 
field entropy, and later is the evolution of the output field statistics. The atomic 
inversion have discussed and plotted against the interaction time Xt. We found that 
it exhibited the conventional Rabi oscillation and collapse-revival for the SDFS. It 
is dependent on the parameters of the used state as squeeze parameter and num- 
ber of photons. We further calculated the field entropy, which is zero for a pure 
state and non-zero for a mixed state. In general , for the atomic radiation system 
the field entropy was found to be non-zero. We have been considered the photon 
number distribution of the output field. We have been obtained the Pegg-Barnett 
phase distribution and also plotted with some parameters. The Q function for some 
parameters has presented analytically and numerically. The Q functin exhibited a 
variety of peaks with corresponding the photon number oscillation. Peak separation 
in the Pegg-Barnett phase distribution and the Q function is associated with the 
onset of the collapse and revivel of the atomic population inversion. The entropy is 
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nearly zero when the Q function has exactly one peak, the greatest separation of Q 
function peaks corresponds to maximum entropy. The effects of squeeze parameter 
and photon number are obvious from all illustrations. The desire to realize physi- 
cally certain specific quantum states such as SDFS's is under current research. It 
is hoped that the SDF states will find applications in the quantum nondemolition 
measurements and quantum optics. They may also find application in experimental 
situations that require low noise sensitivity. 
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